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Spectral approach to the charge carrier transport in p — n-j unctions has aUowed us to revise the 
theoretical results relating to large signal operation of the junction to make them valid for both the 
low and high frequencies ranges. The spectral composition of the external circuit current includes 
both the DC and AC components. The former produces the static current -voltage characteristic 
with the modified Bessel function of zeroth order, Jo(/3V~), while the latter yields the dynamic 
admittance with the modified Bessel function of first order, Ii{l3V^), either of them depend on 
the signal amplitude. Our experimental results are consistent with the theoretical ones if a fitting 
parameter allowing for recombination processes in the depletion layer is taken into account. 
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I. INTRODUCTION 

The conventional theory of p — n-junctions is usually studied in two limiting cases: small-signal/high-frequency 
and large-signal/low-frequency P, 0, H, 0, H, We have found out that it is possible to unify these two cases by 
applying the spectral analysis to calculating the external circuit current for a p—n-diode. General scheme of finding 
the spectral solution to the diffusion equations for the p— n-jimction under arbitrary signal action is developed in 
Sec. [TT] and then applied in Sec. IIIII to calculate the external circuit current. For this work, it is assumed that the 
relaxation times, and the diffusion constants do not change with injection level, or constant effective values are used. 
The nonlinearity considered in this work is related to the doping profile, as it yields a semiconducting p— n-juntion. 
Sections IIVI and IVl deal with the static current -voltage and dynamic characteristics of the diode. First measurements 
to validate the experimental setup is set forth in Sec. IVII Appendix examines the conventional approach (based on 
applying the customary current-voltage characteristic expression to time-varying conditions) to be compared with 
our results obtained below in terms of the spectral approach. 

II. SPECTRAL SOLUTION OF DIFFUSION EQUATIONS 

Theoretical basis of the charge carrier transport in p—n-j unctions is the diffusion-drift equations for minority carriers 
in the n- and p-regions. These equations are derived from the continuity equations together with the current density 
expressions for holes and electrons. Such diffusion-drift equations have the following one-dimensional form 0, 

(a) for holes (with equilibrium and nonequilibrium densities Pn and p{z,t) ) injected into the n-region 

dp djp^ipE) d^jPpp) P-Pn _n, f.-. 

dt ^ dz dz^ Tp ' ^ ' 

(b) for electrons (with equilibrium and nonequilibrium densities Up and n(z,t) ) injected into the p-region 

dn d(n^nE) d^{Dnn) n — r, 



dt dz dz"^ 



0. (2) 



We shall apply Eqs. ([T]) and ^ to study the simplest Shockley's model of the ideal diode [H but operating at 
the nonlinear and multifrequency regime under action of the harmonic signal of arbitrary amplitude. In the case of 
low-level injection in p— n-diodes, one conventionally assumes E = Q for the neutral parts of the p- and n- regions 0, [E] ■ 
Then Eqs. ([T]) and ^ turn into the diffusion equations for the excess concentrations (over their equilibrium values 
Pn and Up) for holes, Ap(z, t) = p(z, t) — pn, and for electrons, An{z, t) = n{z, t) — rip, of the following form: 

d\^^ ^,d'Ap 



l + r,- Ap-L^^.O, (3) 
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(l + -«^jAn-L^,^=0, (4) 

where Lp = ^JD^t^ and L„ = \/DnTn are the diffusion lengths for holes and electrons, respectively. 

For the resonant circuit connection, the total voltage applied to the p— n-diode consists of the DC bias voltage Vq 
and the harmonic signal voltage Vr^ cos cut: 

v{t) = Vo + V^ cos ujt. (5) 

Nonlinearity of electronic processes in the p — n-junction produces frequency harmonics ku so that the required 
solutions of Eqs. ([3]) and (|4|) can be represented in the complex form of Fourier series: 

00 

Ap{z,t)^ ^ Ap,(z)e''="*, (6) 

k — — 00 



An(z,t) = J2 Anfe(z)e^'='^*. (7) 

k — — 00 

Real values of Ap{z, t) and An(z, t) are provided with the following relations for the complex amplitudes: Ap/j — 
Aplj, and Anj, — Anl^,. 

Substitution of the required solutions ([6]) and ([7]) into Eqs. ^ and (jl]) with taking into account the orthogonality 
of harmonics reduces to the following equations for the desired complex amplitudes of harmonics: 

(fApk Apk _ , (fAuk Auk _ 

with 

Lpk = — , and L,, 



y/1 + ikojTp \/l + ikujTn 



General solutions of equations dH) for the p — n-diode with the junction region —Wp < z < Wn have the following 
form: 

(a) for holes injected into the n-region (z > Wn) 

^Pk{z) = C+ exp(^-A,,i^) +C;,exp(^Apfei^'), (9) 



where Apk = apk + ibpk and 



1 / / k(jjT 

apk = ^^l+sJl + ikLUTp)\ bpk^^; (10) 

(b) for electrons injected into the p-region (z < —Wp) 

Anuiz) = C+ exp(^-A„fci±^^ + C*-, exp(^A„fci±^^ , (11) 

where A„fc = a„u + i6„fc and 

1 I / kujTn 
ank = —7= \l I + y I + ykuJTn)'^ , bnk = -Z . (12) 

According to Eqs. (fTO| and (fT2| . always ReAp^ = > and ReA„fe — Unk > 0. So for the p — n-diodes 
with thick base (when dn 3> Lp and dp S> in, where dn and dp are thicknesses of the neutral parts) we 
can take = and C^j, = 0. This eliminates the necessity for boundary conditions on ohmic contacts. Tak- 
ing into account ([9]) and (fTTjl with C^^. — C^j. — 0, the general solutions ([6]) and ([7]) of Eqs. ([3]) and (j4]) assume the form 
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Apiz,t)= J2 Ciexp(-Ap,^-^)e^'=-*, (13) 



k — — oo 



An(z,t)= £ C-,exp('A„fci±i^')e''=-*. (14) 

k — — oo ^ 

These expressions provide for Ap{z,t) — !■ as z — > cx) and An(z,t) ^ as z — > — cx) because of KeApk > and 
ReA„fc > 0. 

The constants and C^^ in Eqs. and (|14p can be found from the conventional injection boundary condi- 
tions 

Ap{Wn,t)^Pnfit) for Z-VK„, (15) 

An(-Wp,i) =nj,/(i) for z ^ -Wp, (16) 
where for the apphed vohage v{t) of the form (5) we have introduced the function 

/(0=exp(^)-l. (17) 

where k = 1.38 x 10~^^ J/K is the Bolztmann constant, and T is the temperature in kelvin. 
Substitution of (fT5|) and into the boundary conditions (|15l) and yields 



^ C+e^'^-*=p„/(0, (18) 



— — OO 



E G«.e^'"*="p/W- (19) 



fc = 



By using the orthogonality property of harmonics in the expansions (|18p and p9p , it is easy to obtain the desired 
constants: 

Cpk=PnFk and C-fe-npFfc, (20) 
where Fk is the Fourier amphtude of the kth harmonic for the function f{t) given by Eq. (jl7p . that is 

Pk = ^ I f{t)e-'''^'diot. (21) 

Substitution of the function (fTT]) into formula (|2T|) gives 

Fo ==/o(/3K.)exp(/3yo)- 1 for fc = 0, (22) 

Fk=F^k=hm^)eMf3Vo) for A: ^ 0. (23) 

Here we have introduced the modified Bessel functions of the first kind of order k {k — 0,±1,±2, . . .) which have 
the following integral representation (see, for example, formula 8.431.5 in [7|): 

7r 

h{PV^) = ^ J e'^^~''°'"^' COS kujtdujt. (24) 



These functions depend on /3V^, where is an amplitude of the signal applied to the n-junction and /3 — q/nT. 
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Consequently, with making allowance for ((20| the general solutions ([13]) and (fT4|) of the diffusion equations ([3]) and 
take the final form of spectral expansions: 



Ap(z,t) = p„ ^^fcexp(-Ap,l-^]e''=-* 



(25) 



(26) 



These expressions allow us to obtain the spectral composition of the current flowing through an external circuit 
connected to the p— n-diode. 



III. EXTERNAL CIRCUIT CURRENT FOR p-n-DIODES 



The initial equation to derive an expression for the diode current is the total current conservation law of the general 
form V • ( jp + j„ + edE/dt) = following from Maxwell's equations. Its one-dimensional form (along the z-direction) 
is 



d_f . dE^ 



0. 



(27) 



The quantity in parentheses of Eq. (|27|) . being independent of z, defines the external circuit current (for the 
cross-section area S of the diode) as a function of time: 



J{t) = jpz{z,t)S + jnz{z,t)S + e^^^^^^S 



(28) 



We have used the usual expressions for the hole and electron current densities [Q] involving the drift and diffusion 
components which underlie the initial diffusion-drift equations ([1]) and All the terms on the right of Eq. ([^5]) 
depend on both z and t, but taken together in any cross section they yield the external circuit current as a function 
of only time. 

Restricting our consideration to the p — n-diodes with low injection (j> <^ n„ but p > Pn so that the excess hole 
density is Ap = p — pn > 0, for n-type, similar conditions hold for p-type), as was noted above, we can assume = 
in neutral parts of the p- and n- regions Q . Then the external circuit current (pS)) is determined only by the diffusion 
currents taken at any cross section, for example, at z = Wn'. 



Jit) 



-qDp 



dpjz, t) 
dz 



S + qD, 



z = W„ 

JpziWn,t)S + Jnz{Wn,t)S. 



dn(z, t) 
dz 



S 



z = W„ 



(29) 



Neglecting recombination processes inside the p — n-junction (which is true if Wn ^ Lp and Wp <C Ln) we can 
write 



jnz{Wn,t) ^ jnz{-Wp,t). 

Substitution of (l30|) into Eq. ([29ll gives the external circuit current: 

Jit) = jpz{Wn,t)S + jn,{-Wp,t)S 



(30) 



dAp{z,t) 
dz 



S + qD, 



z = W„ 



dAn{z,t) 
dz 



S, 



(31) 



-w„ 
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where Ap = p — pn and An = n — Up are the excess concentrations of injected carriers given by Eqs. pK)) and (|26p . 
Inserting these formulas into expression pip , we finally obtain the spectral representation for the external circuit 
current: 

oo oo 

J{t)^Js,p ^feVe'""* + Js.n FkAnkc'"^*. (32) 

k — — oo k — — oo 

Here, the hole and electron contributions into the saturation current of the thick p — n-diode are defined, as it is 
generally accepted [1, Q, in the following form: 

Js,p = — 1^ J and Js,n = — p o- (33) 

Lp Lin 

Expression p2|) contains all the spectral components of the external circuit current, among which the terms with 
numbers fc = (the DC current) and fc = ±1 (the AC current) are of most interest for our subsequent consideration. 



IV. STATIC CURRENT-VOLTAGE CHARACTERISTIC OF p-n-DIODES 

The term in series ((32)) numbered by fc = corresponds to the DC current Jq, which by using expression ([22| for 
-Fb and Apo = A„o = 1 can be written in the form of the static current-voltage characteristic: 

MVo,V^) ^ Js[lo iPV^ ) e''^" - 1 ] . (34) 

Here, lQ{f3V^) is the modified Bessel function of zeroth order depending on the signal amplitude V^. The saturation 
current Jg has the customary form 0, Q 



qPnDp qnpD„ 
b + —— 

Ijp Ijn 



Js — Js.p + Js,n — — f — ~ S H ^ S. (35) 



Expression differs in appearing the Bessel function Io{(3Vr^) from the customary current -voltage characteristic 
given in such well-known books as, for example, [H, Coincidence between them takes place only when Io{pV^) ~ 1, 
which is true for small signals with <ti nT/q. However, an expression similar to our equation (|34p was found 
in a less-known textbook I6| [sl] written in German, whose derivation is adduced for comparison in Appendix (see 
formula fXI)) ). 

The modified Bessel function /o(/3K^) > 1 as a function of the signal amplitude Vr^ mathematically reflects the 
so-called effect of signal rectification. This effect provides a contribution into the DC current Jq from the signal and 
results in up-shifts of the static current -voltage characteristic Jo(Vb) with increasing K.,, as shown in Fig. 1. 



V. DYNAMIC CHARACTERISTICS OF p-n-DIODES 

The first harmonic current contribution to the general expression (I32p for the external current corresponds to terms 
numbered by fc = ±1 and is equal to 

Ji(t) = Jie*"*-H c.c. with Ji =Fi(J,^pApi + J,,„A„i), (36) 

where the quantities Api, A„i, Fi are respectively given by Eqs. pU)) . ((T^), for fc = 1. 

Expressions ([5]) with = 2Vi and pop allow one to introduce for the p— n-diode its dynamic admittance defined 

as 

Y='-^ = Gd + lujCd. (37) 
Vi 

Here, the dynamic conductance Gd and the dynamic {diffusion) capacitance Cd are introduced in a customary 
way 5]. After substituting (|36p into (|37p and some transformations we obtain these quantities as functions of 
frequency: 

Gd(w) =5Go^api(w)^y^ 4- a„i(w)^y^^, (38) 
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FIG. 1: Static current-voltage characteristics for diflFerent values of the signal voltage amplitudes: 1/^=0 (curve 1), fcT/g 
(curve 2), 2KT/q (curve 3), ZnT/q (curve 4), AnT/q (curve 5). 



(39) 



where the quantities flpi, 6pi and a„i, 6„i follow from Eqs. (jlOp and (|12p with A: = 1. 

Formulas (|38p and (|39|) have been derived by inserting the differential conductance introduced from the static 
current-voltage characteristic Jo(Vb, V^) defined as 



Go — 



(40) 



The static current-voltage characteristic, Eq. ([M)) . gives the differential conductance ([3D]) which depends on both 
the bias voltage Vb and the signal amplitude V^: 



r (M 1/ ^ ■^0 + '^'' exp(/3Vb) 
Go(Vo,K.) = — — - — = — /o(/5K.). 



nTjq 



nT/q 



(41) 



Formulas (|38|) and ((39)) contain a combination of the modified Bessel functions /o(/3K^) and /i(/3K.) in the form 
of a factor: 



2 h{pV^) _ gi(y^) 



with 



< 1 



51 (K, 



(42) 



(43) 



These quantities depend on the signal ampHtude K., so that g{Vr^) ~ gi{Vr^) ~ 1 for small signals when ^ 't?'/'? 
and g(V^) ~ 2/(/3K^)— >0 as ^^^oo, whereas here gi{Vr^) is approximated by an exponentially growing function. 
After inserting (gl]) and (gS]) into Eqs. ([38]) and ([391), they take the following form: 



Gd{uj) = giGdo (^api{LLi)'-^ + a„i(w) 



' s.n 



(44) 



\ Qn 



(45) 
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where Q = Qp + Qn with Qp = Js,pTp = qpnLpS and Qn = Js,nTn = qnpLnS. The quantities Gdo and Cdo appearing 
in Eqs. ([M]) and (I45p as functions of Vq are equal to 



qPnl^p Q'f^p^n \ nVn/KT //i'7\ 

CMVo)= 2.r/g =;^(,^- + ^-j^ ■ (^^^ 

Expressions (I46p and (|47p correspond to the diffusion conductance and capacitance at low frequencies (when LUTp,n ^ 
1 and api{uj) — a„i(ti;) ~ 1) for small signals (when Vr^ <C nT/q and gi{Vr^) ~ 1). Our expressions for Gdo and Cdo 
fully coincide with those given in a book ^^i] by Sze (see there formulas (65) and (66) of Chapter 2). 

As seen from PH) and ([iS]) . the frequency dependence of the dynamic conductance Gd{uj) and the diffusion ca- 
pacitance CdiLo) is determined by the functions apilu) and a„i(ti') given by formulas (|10p and (jl2p with fc = 1, 
whereas their dependence on the signal amplitude Vr^ is described by the function gi{V,^) introduced in the form of 
expression (j43| . 

A certain simplification of expressions (|44p and (|45p occurs in the case of the one-sided p'^— n-j unction with highly 
doped emitter when p„ ^ rip, J^^p ^ Js,n, Qp ^ Qn so that Jg ~ Js,p and Q ~ Qp. Then 

gd(^,V~) ^ \/^+^^ + "'''"^' 



CdoiVo) 

where from Eqs. (^5)1 and ([i7|) it follows that 

GdoiVo) 

Frequency dependencies given by expressions (|48p and are plotted in Fig. 2. The plots differ from the similar 
curves shown by Sze Q (see there Fig. 23 of Chapter 2) in that they take into account the signal amplitude influence 
owing to the factor gi{Vr^) = Ii{f3Vr^)/{(3V^/2). For small signals (when f3Vr^ <C 1 and (7i(K.) ~ 1) our expressions 
and ([^5]) assume the form obtained by Sze [5] and are depicted by curves 1 in Fig. 2. The low-frequency values 
(when ojTp <C 1) of the quantities (|48| and ((49|) for arbitrary signals are equal to 

GdiO,V^) ^ Cd{0,V^) ^ hiPV^) ^ , . 
GdoiVo) Cdo(^o) /3^~/2 "^'^ 

From here it follows that the low-frequency value (denoted by zero in parentheses that means = 0) of the dynamic 
conductance is described by the following expression: 

Gd{Q.V^) = Gd^{Vo)gi{V^). (51) 

Eq. (|5ip coincides with formula (|A5p obtained from the conventional approach based on using the customary 
current-voltage characteristic for time-varying conditions. 

The conventional approach is not able to provide the frequency dependence not only for Gd{io) but also for Gd{oj) 
such as given by Eq. and which are valid for signals of arbitrary amplitude. The similar frequency dependency 
known from the published literature was obtained only for small signals when Vn^<^KT / q. Such a result can be found, 
for example, in Sec. 2.4.4 of a book by Sze [5], where the curves Gd{oj)/Gdo and Cd{uj)/Cdo given in Fig. 23 correspond 
to our curves 1 in Fig. 2 calculated for Vr^ = 0. Other our curves denoted by numbers 2, 3, 4, 5 in Fig. 2 are novel 
since they demonstrate a dependence of the dynamic conductance and capacitance on the AC signal amplitude V^, 
which was previously missed. It is this dependence that is of interest to check experimentally. 



V2 Ii{f3V^) 
1 + JI+uj^'t^ 



2 ^ . JseMPVo) 



(49) 



(50) 
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FIG. 2: Frequency dependencies of the normalized dynamic conductance Gd{i^) /Gdo and diffusion capacitance Cd{i^) /Cdo for 
different values of the signal voltage amplitudes: = (curve 1), nT/q (curve 2), 2nT/q (curve 3), SuT/q (curve 4), AkT/q 
(curve 5). 



VI. EXPERIMENTAL VALIDATION 



The above theoretical results have been derived from the ideal Shockley's model which have a few simplifying 
assumptions. Among them the most essential ones are: (i) the large thickness of the neutral regions as compared 
with the minority-carrier diffusion lengths, i. e., d„ 3> Lp and dp ^ L„; (ii) the neglect of recombination effects in the 
depletion layers owing to inequalities Wn ^ Lp and Wp ^ L„. Yet, such assumptions are rarely verified in practical 
conditions and should be taken into account, as applied to real devices used in our experiments. 

The practically used p'^— n-diode with the one-sided injection has the thin n-base of thickness (i„ compared to Lp. 
In this case, it is necessary to apply the boundary condition taking into account surface recombination on a metallic 
contact of the base located at z = Wn + dn (see Eq. (108) of Chapter 1 in Ref. [5j]): 



dAp{z, t) 



dz 



= -|^Ap(l^„+d„,t), 



where Sp is the speed of surface recombination of holes (for ohmic contacts Sp oo). The additional boundary 
condition changes only the saturation current by the factor Ap-. 

_qPnDp_ . . jspLp/Dp) + tanh((14^„ + d„)/Lp) 

'^'^ Lp ^'^"^ "^''^ \ + {spLp/Dp)ts.nH{Wn + dn)lLp)- ^^'^ 

As follows from Eq. ((50|) . the measurement of the small-signal low- frequency conductance Gdo — PJs exp(/3Vo) will 
give us the saturation current (|52p with allowance for the surface recombination on metallic contact of the base. 

To take into account the bulk recombination processes in the depletion layer, let us follow a phenomenological ap- 
proach suggested by Sze |5[ and introduce the empirical nonideality factor n so as to provide the following replacement: 

/5 ^ ^ ^ ^ ^ /?„. (53) 

Values of the factor n lie between 1 and 2 [5]: (a) n = 1 when contribution of the bulk recombination processes is 
negligibly small, (b) n ^ 2 when the recombination current dominates over the diffusion one. 

By using the replacement ()53|) . expression ()5ip for the low- frequency dynamic conductance can be rewritten in the 
following corrected form: 

Gd{0,V^) = GdoiVo)gn{V^), (54) 
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FIG. 3: Experimental setup with the measured diode D and reference capacitor C. 



where the correcting function gniV^) is defined as 



(55) 



In general, the modified quantity /3„ in expression ()55p takes into account not only a contribution from the bulk 
recombination processes by means of the nonideality factor n but also that from a priori unknown temperature T of 
the p— n-junction under experimental investigation. Hence, the product nT can be used as a fitting parameter to 
adjust the theoretical relations ((54|) - ([55)) with experimental results obtained below. 

To verify our theory, we built a simple apparatus using the dual-phase DSP lock-in amplifier, Stanford Research 
Systems model SR830. For the low frequency measurement, the low noise electrometer grade operational amplifier, 
Burr-Brown OPA128JM, is used. For the high-frequency measurements, the same lock-in amplifier, combined with 
a mixer, and high-frequency op-amps can be used, see Fig. 3. The AC voltage is applied to the diode with no bias 
(Vb = 0). The lock-in amplifier has two displays — X and Y, which give the root-mean-square (rms) value of output 
signal at the excitation frequency uu. It is easy to show that the lock- in output in the X and Y displays is equal to 



X 



-RV„ 



Gd{Lo) cos(0 - (t)o) + ujCd{Lo) sin((?!) - 0o) 



(56) 



Gdiuj) sm{(j) - (j)Q) + uCdiuj) cos(0 - (/iq) 



(57) 



Here i? is a feedback resistor, Vm is the rms value of the lock- in oscillator voltage {Vin = is an internal 

phase of the lock-in oscillator signal connected to the phase detector, and 0o is a phase of the signal after passing 
through an external circuit. From (|56p and ([57]) it follows that the two displays are 7r/2 out of phase from each other. 

To get correct experimental values, it is necessary to adjust the phase (j) of the lock-in local oscillator so that the 
X-display would be used for the conductance voltage \V^i\ = RVmGd and the Y-display for the capacitance voltage 
iKmtl — R^in^Cd- As follows from expressions and ([57)) . it can be realized only if = 0o- The phase (/) is adjusted 
with the reference capacitor C which is placed between the local oscillator and the minus input of the operational 
amplifier instead of diode as shown in Fig. 3. The phase adjustment is carried out until zero voltage is observed 
in the conductance X-display, i. e., when (f) = (f)Q. Then Eqs. ((56)) and ([57]) yield the required results: 



Gd{iJ,V^) - - 



1 \Kit 



R V„. 



1 \V^ I 



R V„ 



(58) 



After the adjustment, the reference capacitor C is replaced by the measured diode D and the measurements are 
performed by varying the applied AC voltage Vr^ with Vq = 0. The operating frequency of the lock-in internal 
oscillator was chosen IkHz to surely provide the relation uTp <C 1 underlying the initial theoretical expression (I5ip . 
To extract the conductance, one must normalize the data. 
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FIG. 4: Theoretical and experimental dependencies of the normalized dynamic conductance Gd{0, V^)/Gdo on the peak voltage 
of AC signal for the 1N914B diode with no DC bias (Vo = 0). 

For our experiments, we have employed a common high frequency diode 1N914B. The experimental results for 
the above-mentioned diode are presented in Fig. 4 which demonstrates a quite satisfactory agreement between the 
experimental data and our theoretical results, in the low frequency range. All resistors and capacitors employed in the 
experiment were measured with the Stanford Research Systems model SR720 LCR meter. The measurement error is 
found to be below 1%. In an additional paper, the detailed experimental results will be presented. 

As seen from Eqs. (j54p and (|55p . the small voltage measurement (when gn= 1) gives a value of the low-frequency 
dynamic conductance Gdo defined by formula ((50)) . This value has proved to be equal Gdo — 6.4 x 10~^mho. From 
here it follows that Jg — Gdo/Pn ~ 31 nA, which is in agreement with the reverse current for the 1N914B diode given 
by manufacturers. 

From the experimental curve plotted in Fig. 4 it follows that the modified factor /3„ = q/nuT should be equal to 
20.6 for the 1N914B diode in order to fit the theoretical expression (|54p . Hence, the fitting parameter nT is equal 
to 560 K, which for the operating temperature T — 400 K provides the nonideality factor n — 1.4. Therefore, the 
measured diode operates in a regime when the diffusion current slightly dominates over the recombination current. 

VII. CONCLUSION 

Spectral approach to the theory of p—n- junctions has allowed us to take into account the effect of large signal at 
both the low and high frequencies as compared to t~^. This approach is based on the known diffusion equations for 
injected minority carriers, which is a standard practice for semiconductor electronics. 

The only specific feature distinguishing our approach from the conventional one given for comparison in Appendix 
is related to the initial representation of the desired carrier concentrations in the form of Fourier expansion over 
frequency harmonics. Such harmonics are produced by nonlinear processes in the p— n-junction when the sufhcicntly 
large AC voltage Vi — is applied to the junction together with the DC bias voltage Vq. As a result, the spectra of 

both the excess concentration of injected carriers and the external circuit current have been derived (see formulas (j25p . 
p6|. and ((32)) ). The use of appropriate terms in Fourier series has given rise to expressions for the DC component Jq 
and AC component Ji of the external circuit current. The former determine the static current- voltage characteristic 
Jo(Vb) (formula ((M| ) and from the latter follows the dynamic admittance Y{uj) = Ji/Vi (formulas (|37 p - ((55)) ). These 
expressions have proved to be dependent not only on Vq but also on Vi , in the case of the nonlinear (in signal) regime 
of operation of the p— n-junction. 

The results are corroborated by experimental verification, as follows from Fig. 4. Detailed experimental results will 
be presented in a separated paper. 
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APPENDIX A: CONVENTIONAL APPROACH TO ANALYSIS OF THE CURRENT-VOLTAGE 

CHARACTERISTIC FOR p n-DIODES 



Conventional theory of the p— n-diode operation under the AC signal of arbitrary amplitude is based on generalizing 
the known expression for the static current-voltage characteristic to time- varying conditions: 



(Al) 



Here the total voltage v{t) applied to the diode includes the DC bias voltage Vq and the harmonic signal voltage 
Vr^ cos (see Eq. ([5])). After substituting ([5]) into Eq. (jAip . it is necessary to use the following Fourier series (see 
also formula 8.511.4 in 0): 



f3Vr^ COS LUt 



lo (/3K. ) + 2/i {PV^ )cosLut + 2l2il3V^) cos 2ujt + 



(A2) 



where /o,/i,/2, . 
sion (|Aip yields 



are the modified Bessel functions of first kind. Insertion of the Fourier series (IA2p into expres- 



At) 



'2l2{(3V^)e^^° cos2ujt 



(A3) 



Expression (|A3p is the Fourier expansion of the diode current whose DC component gives the static current -voltage 
characteristic 



Jo(Vo,y.) = J. /o(/31/^)e^^" -1 



(A4) 



which is valid for the harmonic signal of arbitrary amplitude . The first harmonic of the current (jASp allows one 
to obtain the dynamic conductance 



GdiVo,V^) 



KT/q l3V^/2 



Gdo{Vo)gi{Vr.), 



(A5) 



where GdoiVo) and gi{Vr^) are given by formulas (|46| and (l43l) . 

Expression (jA5p resulting from the conventional approach describes solely the low-frequency dynamic conductance 
Gd (cf. Eq. (j5ip ). whose frequency dependence can be obtained only by using the spectral approach, as it is done in 
Sec. |Vl Moreover, our spectral approach yields not only Gd{i^) but also Cd(ijj) for arbitrary signals. In the case of 
small signals, the similar frequency dependencies are derived in the literature (for example, see HI), but they never 
follow from the current -voltage characteristic in the form of Eq. (|Aip taken as an initial point for derivation. 
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